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Abstract: We present some recent developments in the theory of particle acceleration at shock fronts in 
the presence of dynamical reaction of the accelerated particles and self-generation of magnetic field due 
to streaming instability. The spectra of accelerated particles, the velocity, magnetic field and temperature 
profiles can be calculated in this approach anywhere in the precursor and in the downstream region. The 
implications for the origin of cosmic rays and for the phenomenology of supernova remnants will be 
discussed. 
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Introduction 

The paradigm of supernova remnants as sources of 
the galactic component of cosmic rays heavily re- 
lies on the fact that the mechanism of particle ac- 
celeration at collisionless shocks should take place 
in the nonlinear regime. The nonlinearity mani- 
fests itself in at least two respects: 1) the efficiency 
of acceleration necessary to explain the origin of 
cosmic rays needs to be large enough that the dy- 
namical reaction of the accelerated particles is not 
negligible; 2) the accelerated particles are respon- 
sible for the self-generation of the magnetic distur- 
bances which in turn scatter the particles thereby 
allowing their acceleration to the maximum ener- 
gies observed in cosmic ray data, most notably the 
KASCADE data in the knee region [1]. 

A special emphasis should be given to the fact that 
protons may be energized to energies of 10^ — 10^ 
GeV only if the magnetic field in the shock vicinity 
is amplified by a factor few hundreds and reorga- 
nized in the form of a flat power spectrum which 
may lead to Bohm diffusion. Such amplification 
may take place through streaming instability in- 
duced by the super-Alfvenic drift of accelerated 
particles in the frame of the plasma upstream of 
the shock [2, 3, 4]. Magnetic field amplification 
may also take place due to firehose instability [5]. 



A comprehensive theory of particle acceleration in 
the nonlinear regime, with both dynamical reaction 
of the accelerated particles and magnetic amplifi- 
cation taken into account has been recently formu- 
lated in [6, 7]. In [8] the calculation of the maxi- 
mum momentum for the nonlinear regime was pre- 
sented. 

Here we illustrate the basic formalism and the phe- 
nomenological implications of the theoretical ap- 
proach of [7]. 

The formalism 

The two basic equations needed in this section are 
the equation of conservation of momentum and the 
transport equation for the accelerated particles. In 
the upstream plasma, conservation of momentum 
reads: 



Ux) = 1+ 
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where ^c{x) = PcB.{x)/poul and U{x) = 
u{x) I uq and we used conservation of mass pqUq = 
p{x)u{x) (here and uq refer to the density and 
plasma velocity at upstream infinity, while p{x) 
and u{x) are the density and velocity at the loca- 
tion X upstream. A/q is the sonic Mach number 
at upstream infinity). The pressure in the form of 
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accelerated particles is defined as 

Pcr{x) = \ T"""^ dp 4np^vip)f{x,p), (2) 

J Pin, 

and f(x,p) is the distribution function of acceler- 
ated particles. Herepinj andpmaa; are the injection 
and maximum momentum. The function / van- 
ishes at upstream infinity. The distribution function 
satisfies the following transport equation in the ref- 
erence frame of the shock: 
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[6] and [9] showed that an excellent approximation 
to the solution f{x,p) has the form 



f{x,p) = /o(p)exp 
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where fa{p) = f{x = 0,p) is the cosmic rays' 
distribution function at the shock and q{p) = 
— is its local slope in momentum space. 

The function /o {p) can be written in a very general 
way as found by [10]: 




exp 



Here we introduced the function Up{p) = Up/uo, 
with 



Up = Ul 



1 /-o 

fo{p) i-oc 



dx{du/dx)f{x,p) , (4) 



where ui is the fluid velocity immedi- 
ately upstream (at x = 0~). We used 

Qix,P) = '-^^B^Sip - p.nMx), With 

ngas,i = noRtot/Rsub the gas density imme- 
diately upstream (x = 0^) and 77 the fraction 
of the particles crossing the shock which are 
going to take part in the acceleration process. 
In the expressions above we also introduced the 
compression factor at the subshock Rgub = u\/u2 
and the total compression factor Rtot = uo/u2. 



If the upstream plasma only evolves adiabatically, 
the two compression factors are related through 
the following expression ([10]): 



Rt, 
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where Mq is the Mach number of the fluid at up- 
stream infinity and 7^ is the ratio of specific heats 
for the fluid. The parameter 77 in Eq. 4 contains the 
very important information about the injection of 
particles from the thermal pool. The injection is 
modelled as proposed in [1 1]: 
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(6) 



Here ^ is a parameter that identifies the injec- 
tion momentum as a multiple of the momentum 
of the thermal particles in the downstream sec- 
tion (pinj = £,Pth,2)- The latter is calculated self- 
consistently from the Rankine-Hugoniot relations 
at the subshock. For the numerical calculations 
that follow we always use ^ = 3.5, that corre- 
sponds to a fraction of order 10~^ of the particles 
crossing the shock to be injected in the accelerator 

The scattering properties of the background plasma 
are described by the scalar function D{p), the dif- 
fusion coefficient. Once J-{x, k) is known, the dif- 
fusion coefficient is known in turn ([2]): 



D{x,p) = - 



(7) 



From the latter equation, where vl stands for the 
Larmor radius of particles of momentum p, it is 
clear that the diffusion coefficient tends to Bohm's 
expression for ^ — > 1. The expected saturation 
level for the overall energy density of the perturbed 
magnetic field can be easily evaluated from the fact 
that 
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J^(k, x) = VA 



dPcR 
dx 



(8) 



Integration of this equation is straightforward 
when non-linear effects on the fluid are neglected 
so that u and va are both spatially constant. One 
obtains SB'^/Stt =- {va/u)Pcr, or, in terms of 
amplification of the ambient magnetic field: 
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(9) 
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Figure 1: Spectrum and slope at the shock loca- 
tion as functions of energy forpmox = lO^mc and 
magnetic field at upstream infinity Bo = l/iC The 
curves refer to Mach numbers at upstream infin- 
ity ranging from Mq = 4 to Mq = 200: dotted 
for Mo = 4, dashed for Mo = 10, dot-dashed for 
Mq = 50, solid for Mq = 100 and dot-dot-dashed 
for Mo = 200. 

with Ma = uq/va the Alfvenic Mach number. 

It is worth stressing that for Pcr/pqUq 1 and 
Ma ^ 1, the predicted amplification of the mag- 
netic field exceeds unity. CLearly this implies that 
the quasi-Unear theory used here loses vahdity and 
that a more accurate description, though very diffi- 
cult to achieve should be sought. 

The set of equations of conservation of momentum 
and transport equation with a diffusion coefficient 
determined as described above can be solved by 
using the iterative procedure described in detail in 
[6, 7]. 

Results 

The spectra of the accelerated particles for Mach 
numbers at upstream infinity ranging from Mq = 4 
to Mq = 200 are shown in Fig. 1 for a background 
magnetic field at upstream infinity Bq = IfiG (the 
result is however independent of the strength of the 
background magnetic field). In the bottom part of 
the same figure we plot the slope of the spectrum 
as a function of momentum. 

For low Mach numbers and at given Pmax the mod- 
ification of the shock due to the reaction of the 



accelerated particles is small (see for instance the 
case Mo = 4). For the strongly modified case (e.g. 
Mq = 200) the asymptotic spectrum of the accel- 
erated particles is very flat, tending to with 
a = 3.1 - 3.2 forp Pmax- 

We need to comment on the issue of the shape 
of the spectrum of cosmic rays accelerated in su- 
pernova remnants: the spectra illustrated in Fig. 
1 are the spectra in proximity of the shock sur- 
face and therefore the ones which are relevant for 
the calculation of the spectra of secondary radia- 
tion produced by the accelerated particles. In gen- 
eral the spectra that are observed by an observer 
far upstream are more complex to determine: at 
each time during the supernova evolution parti- 
cles at the maximum momentum can escape to up- 
stream infinity (this is a peculiar aspect of non- 
Unear theory) carrying away a sizeable fraction 
of the total energy (due to the flat spectra). At 
each time the instantaneous spectrum that escapes 
to upstream infinity is a narrow function centered 
around Pmax(t), where t is the age of the rem- 
nant. If the maximum momentum decreases with 
time (in the Sedov phase this is the case) then a 
spectrum is built due to the overlap of many delta- 
function-like spectra leaving the acceleration re- 
gion from upstream. In the simple model con- 
sidered in [12] this overlap leads to a power law 
spectrum despite the fact that the spectrum at 
the shock may be concave as illustrated in Fig. 1. 
In addition to these particles that leave the system 
from upstream, the distant observer will also mea- 
sure the cosmic ray spectrum which is kept in the 
supernova shell and is eventually liberated at later 
times, possibly suffering adiabatic energy losses. 
The spectrum observed at the Earth is likely to be 
a complex superposition of the particles escaping 
from upstream, those leaving the system at the end 
of the supernova evolution, and summed over all 
supernova events occurred during a confinement 
time of cosmic rays in the Galaxy. 

The diffusion coefficient associated with the self- 
generated waves is given by Eq. 7. We plot this 
diffusion coefficient at the shock location in Fig. 2 
for Mach numbers Mq = 10 (dashed Unes) and 
Mo = 100 (solid lines). For comparison, we also 
plot the corresponding Bohm diffusion coefficient 
Db{p) oc v{p)p in the unperturbed magnetic field 
Bo, for Bo = l/iG and Bq = lO/xG. 
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Figure 2: The self-generated diffusion coefficient 
at the shock location a; = 0~ as a function of the 
particle momentum for Mach numbers Mq = 10 
(dotted fine), Mq = 100 (dashed line) and Mq = 
200 (solid Une). Also plotted is the Bohm diffusion 
coefficient corresponding to Bq = l^G (solid line 
with triangles) and Bq = lO/xG (solid line with 
diamonds). The y-axis is in units of cm^s~^. 

As stressed above, in the regime we considered, the 
fluctuations in the magnetic field become strongly 
non linear. The dynamical role of the amplified 
field remains however negligible as the highest val- 
ues of 6B^ /SttpqUq, reached close to the shock 
front, are of the order of 10~^ - 10~^. 

Conclusions 

We developed a mathematical formalism that al- 
lows us to calculate the spectrum of particles ac- 
celerated at a collisionless non-relativistic shock 
taking into acount both the nonlinear dynamical 
reaction of the accelerated particles and the mag- 
netic field amplification by streaming instability. 
The approach has also recently been generalized 
to allow us to determine self-consistently the max- 
imum momentum reached by the particles in this 
fully nonlinear acceleration regime [8]. 
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